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Abstract 

In this note we present a mathematical framework for a rigorous approach to a 
common track fit for trackers located in the inner region of the ZEUS detector. The 
approach makes use of the Kalman filter and offers a rigorous treatment of magnetic 
field inhomogeneity, multiple scattering and energy loss. We describe mathemat- 
ical details of the implementation of the Kalman filter technique with a reduced 
amount of computations for a cylindrical drift chamber, barrel and forward silicon 
strip detectors and a forward straw drift chamber. Options with homogeneous and 
inhomogeneous field are discussed. The fitting of tracks in one ZEUS event takes 
about of 20 ms on standard PC. 
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1 Introduction 



The ZEUS experiment pQ was operated at the electron-proton colhder HERA at DESY 
until 2007. The ZEUS detector was a sophisticated, multi-component tool for studying 
particle reactions provided by electron-proton collisions with an energy 27.5 GeV and 
920 GeV,respectively. The inner tracking components of the ZEUS detector were: the 
silicon strip Micro Vertex Detector [2] with barrel (BMVD) and forward (FMVD) parts; 
the Central Tracking Detector (CTD) [3] consisting of the cylindrical drift chamber; the 
Forward Tracking Device (FTD) [T] and the forward Straw- Tube Tracker (STT) |4j. The 
MVD was located in the vicinity of interaction point, inside of the CTD. 

The magnetic field in the central region of the ZEUS detector was produced by a thin 
superconducting solenoid. The field had a strength of 14.3 kCauss at the center and was 
directed parallel to the proton beam. The barrel MVD and CTD were located in the 
field which was almost homogeneous with a small radial component far from the center. 
Forward trackers were placed outside of the solenoid or close to its edge where the field 
is inhomogeneous. 

We consider a mathematical framework for a rigorous approach to a common track 
fit, which can be performed with tracks including all inner tracking components or with 
any combination of them. The approach offers a rigorous treatment of field inhomo- 
geneity, multiple scattering and energy loss. The track fitting procedure makes use of 
the Kalman filter technique and we discuss how to optimize computations and make the 
fitting procedure fast. 

2 Overview of the tracker layout 

The ZEUS coordinate system is a right-handed Cartesian system, with the z-axis pointing 
in the proton beam direction (forward) and the x-axis pointing to the center of the HERA 
ring. The coordinate origin is at the nominal interaction point. 

The barrel (BMVD) and forward (FMVD) section of the MVD includes 600 and 112 
sensors, respectively [2]. A sensor is a silicon single-sided strip detector with a readout 
pitch of 120 fim which includes five innermost strips for capacitive charge division. The 
ZEUS MVD has 307,200 and 53,730 readout channels in the barrel and forward sections, 
respectively. 

The barrel section, centered at the interaction point, is about 63 cm long. The silicon 
sensors are arranged in three concentric cylindrical layers with radii about 5 cm, 8 cm and 
12 cm. Two back to back sensors in a layer provide measurements of nominal r — (j) and 
z position. The FMVD is composed of four transverse disks of 14 wedges each, which 
extend the angular coverage down to 7° from the beam line. Each wedge has two sensor 
layers separated by approximately 8 mm in 2;-direction. They are mounted back to back, 
such that the angle between strips is 2 x 13°. 

The CTD |2] is a cylindrical drift chamber, with a sensitive volume approximately 2m 
in length and 0.4 (1.6m) in inner (outer) diameter. The CTD wires are arranged into nine 
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concentric super layers numbered consecutively from the inside out. The odd- numbered 
superlayers have sense wires running parallel to the chamber axis (i.e. z-axis) while 
those in the even-numbered superlayers have a 5° stereo angle. We denote sense wires 
in corresponding superlayers as axial and stereo, respectively. Each super layer contains 
eight sense wire layers - there are 4608 sense wires in total. A set of eight sense wires is 
surrounded by field wires, azimuthally dividing a superlayer into cells of polygonal shape. 
Each sense wire is read out by a flash ADC and, finally a drift distance is evaluated for a 
hit. All axial wires in superlayer one and the odd numbered wires in superlayer three and 
five (in total 704 wires) are additionally equipped with the z-by-timing system, which 
measures z position of a hit. 

The STT uses straw drift chambers with 7.5 mm diameter capton tubes of varying 
length from 20 cm to 75 cm. There are in total 10,944 wires in 48 wedge shaped sectors. 
Each wedge covers an azimuthal span of 60°. Each sector consists of 3 layers of straws 
perpendicular to the 2;-axis. A track crossing the STT nominally dehvers 24 drift time 
measurements. 



3 Track Models and Likelihood Functions in a Multi- 
Component Tracker 



The hkelihood function of a track measurement has a meaning regardless of the details 
of any fitting method. The maximum-likelihood estimator is efficient in the sense that 
no other unbiased estimator has smaller variances. A track model which is appropriate 
for the likelihood function, together with a given method of track fit, may produce an 
efficient estimate of parameters. A general point of view of the information delivered by a 
tracker can help to interpret behavior of variances of fitted parameters and hit residuals. 

We can model a multi-component tracker by a set of track detecting elements and 
intermediate blocks of passive material, which are located in a static magnetic field. 
Track parameters in the detector element k are described by a vector Xk- For the case of 
a three-dimensional fit, the dimension of the vector, x^, is 5. The track measurement in 
the tracker element k, i.e. the k^^ hit, is a vector denoted by ruk. In general is the 
vector with its dimension corresponding to that of the tracking element. For example, 
rrik has only 1 coordinate for a silicon strip of the MVD, a drift tube of the STT or 
a stereo wire of the CTD and 2 coordinates (drift time and z position of a hit) for an 
axial wire of the CTD which is additionally equipped with the z-by-timing system. The 
measurement error can be described by the covariance matrix Vk- We approximate the 
probability (density) of the measurement given the vector of track parameters Xk 




(1) 



by a Gaussian function with the mean value (mfc) and covariance matrix V^: 



G{mk\{mk);Vk) = C(Vfe)exp{ - -(m^ - {mk)fVf. ^{ruk - (mfe))}, (2) 
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where C{Vk) is a normalization constant. An operator Hk projects the actual vector Xk 
into the space of measurement: 

{mk) = HkXk. (3) 

Suppose that we are interested in the track parameters at the beginning of track, xi. 
The likelihood function takes the form of a product: 

N 

L{mi,m2,...,mN\xi)=P{x2,...,XN\xi)- Y[G{mk\HkXk;Vk). (4) 

k=l 

The first term is the probability for a particle to pass through the points X2, ■■.,X]y given 
the parameters Xi at the beginning and the second one is the probability to obtain the 
measurements, mi, ...,mjv, while measuring the points in the space of track parameters 
xi, ...,xn of the real (not the mean) trajectory. The probability, P{x2, ...,xn\xi), can be 
approximated by a Gaussian distribution 

P(X2, ...,Xn\Xi) = G{X2, ...,XN\{x2iXi)), {xNiXi))] T,{xi)). (5) 

The mean trajectory is defined as: 

(xfc(xi)) = J^kXi, (6) 

where the operator J-'k swims track parameters xi into the detector element k. The track 
model may be described as a continuous curve for the mean trajectory with fluctuations 
of actual parameters Xk with respect to the mean trajectory, 

Vk{xi) = Xk- J^kXi- (7) 

The fluctuation, T>k{xi), accumulates the effect of multiple scattering on the pass from 
the beginning of the track to the given element. Vectors {Vk{xi)} are correlated and, 
therefore, matrix has dense structure (many non-zero elements). We can combine 

Gaussian functions from (j4]) and ([5]): 

L{mi,m2, ...,m7v|xi) = G'(mi,m2, ...,mN\HiXi, H2J^2Xi, HnJ^nXi, M{xi)), (8) 

where the non-diagonal covariance matrix Ai{xi) has dimension equal to the sum of 
dimensions of all measurements {rrik}. The dimension of the Ad may be of order 10^ for 
modern tracking detectors. Maximization of the likelihood function of Gaussian type, 
i.e. least square fitting with large non-diagonal covariance matrix Ai, requires a lot 
of computations, although not more than 5 parameters are fitted. Because of large 
computing time, the model is not convenient for a track fitting in a multi-component 
tracker. But the model includes a small number of fitted parameters, and is suitable 
for a subsequent update of detector alignment parameters [5] , where an expansion of hit 
residuals w.r.t. fitted parameters is needed. 

A charged particle traversing a medium can be described by a stochastic process with 
the Markov property and, therefore, the conditional probability distribution of future 
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states depends only upon the present state and not on any past states. The probabihty 
function for a particle to pass through the points x2, ...^xn in ([5]) can be rewritten as: 

N N 

P{X2, Xn\Xi) = Yl PiXk\Xk-l) = n G{Xk\{Xk{Xk-l))]Qk{Xk-l))- (9) 
k=2 k=2 

We approximate the conditional probability (density), P{xk\xk-i), for track parame- 
ters Xk, given the parameters in the previous state Xk-i, by the Gaussian distribution, 
G{xk\{xk{xk-i)); Qk{xk-i)) with the mean {xk{xk-i)) and covariance matrix Qk{xk-i)- 
The mean trajectory in the tracking element k is 

{xk{xk-i)) = FkXk-i- (10) 

The operator swims track parameters x^-i into the detector element k according to 
the equations of motion. 

Suppose that we are interested in track parameters in all points of track measurement, 
i.e. Xi,X2, ...,xn- The likelihood function takes a form: 

N 

L(mi, ...,mN\xi, ...,xn) = G{mi\HiXi; Vi) ■ G{xk\FkXk-i; Qk) G{mk\HkXk] V^), 

k=2 

(11) 

with Gaussian functions 

G{nik\HkXk] Vk) = C(Vfc)exp{ - ^(m^ - HkXkfVj;^{nik - HkXk)] (12) 

and 

G{xk\FkXk-.i] Qk) = C{Qk) exp | - ^{xk - FkXk-if Ql^{xk - FkXk^i)], (13) 

where CiVk) and C{Qk) are normalization constants. 

The model for the total track is not a continuous curve, but consists of — 1 contin- 
uous segments. A variation of track parameters in the point of discontinuity 

5k = Xk- FkXk-i (14) 

describes the effect of multiple scattering on the pass from the the previous element k — 1 
to the element k. Vectors {6k} are uncorrelated. A spread of the 6k is defined by the 
covariance matrix Qk- 

The maximum likelihood estimation of parameters {xk} satisfies the system of equa- 
tions 

If operators Fk and Hk are non-linear (e.g. in magnetic field) then the latter equa- 
tions are non-linear also. The problem can be solved iteratively using the well known 
method of linearization of operations ([3]) and (ITOl) . Anyhow we can regard the functional. 
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d{—lnL)/dxJ, as a linear form w.r.t. vectors of estimated parameters {xk}- The vector, 
Xk, associates in ffTTl) only with vectors in neighboring data points k — 1 and k + 1 and, 
therefore, the linear form d{—\nL)/dx1 includes only 3 terms with vectors Xk-i,Xk and 
3^fc+i05 respectively. Finally, the system IHM looks as 
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where submatrices related to points i,j, 

^ _ d\-\nL) 

dxfdxj 

are parts of the information matrix. The sparse (with many zero elements), band struc- 
ture of the information matrix can be exploited to reduce computations drastically. This 
can be achieved by using either a dedicated algorithm of matrix inversion [6J, or else (e.g. 
in the broken lines fit [7]) by the matrix (Cholesky) decomposition into a unit triangle, 
U, and a diagonal, D, matrix 

U D f/^x = r 
which requires two steps to solve for x: 

Uy = r and D U^x = y . 



The track model based on relations ( jTOl - UM is well suited also for an implementation 
of the progressive track fit by the method [S| or for the application of the Kalman filter 
formalism [9]. Both methods are rather economical regarding computing time because 
they include operations with matrices of maximal size 5 by 5 for each hit. 



4 Application of the Kalman filter technique to track 
fitting 

In [2] it was shown that an appropriate mathematical framework for the iterative proce- 
dure of track fitting is the theory of linear filtering, in particular the Kalman filter [10] . 
To consider the mathematical framework of a Kalman filter, we try to follow the notation 
used in [TT]. In the following we describe a case with a linear system and a non-linear 
system will be discussed at Subsec. 14.21 

^Differentiating the latter linear form with respect to Xi isolates a coefficient in a corresponding linear 
term. 
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4.1 Linear Model 



The Kalman filter proceeds progressively from one measurement to the next and improves 
the knowledge about the particle trajectory by updating the track parameters with each 
new measurement. The system state vector (track parameters) after inclusion of k — 1 
measurements is denoted by Xk-i, and its covariance matrix by Cfc-i- The state vector 
and its covariance matrix are propagated to the location of the next measurement with 
the prediction equations: 

x^-' = FkXk-i, (17) 

and 

^ = FkCk-iF^ + Qk, (18) 

where Fk is the transport matrix and Qk denotes the covariance matrix of the process 
noise, which occurs due to the random perturbation of the particle's trajectory. 

The measurement of the vector x'^"^ and its covariance matrix are denoted by and 
Vfe, respectively. The expectec? measurement is described by the projection matrix H^. 
The estimated residuals are 

rl-'^rrik-Hkxl-' (19) 
and its covariance matrix become: 

Rt'^Vk + HkCfHl (20) 

The updating of the system state vector after inclusion of the measurement k is defined 
by the filter equations: 

Xk = xl~' + Kkrt\ (21) 
Cfe = (1 - Kk Hk) Cfe~\ 
with the filtered residuals and its covariance matrix 



Tk 



[1 - Hk Kk) rt\ Rk = {l- Hk Kk) Vk = Vk- HkCkH^ (22) 



The matrix, K^, is called the filtering (gain ) matrix. The filtered state vector is pulled 
towards the measurement and, therefore the quadratic mean of the filtered residual is 
smaller than the measurement error. The increment after the filtering of the state 
vector is given by: 

xl = rlRk^Tk. 

The track parameters after the filtering procedure are known with optimal precision 
only at the last point of the fit. The smoothing part of the Kalman filter is a very useful 
complement, which solves the problem of optimal parameter estimation at every point of 
the trajectory. The smoothing is also a recursive procedure which proceeds step by step 
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in the direction opposite to that of the filter with the smoother equations: 



— Ck F'^+i (C'fc+i) ^-i 
C]: = C, + A,{C^^,-Ct^,)Al, (23) 

Rl = Rk- HkAk{Ck_^_i - Ck^^)AlHl = Vk- HkClHl. 

The smoothed state vector, x^, is more precise, because it includes information from 
all measurements. The variance of the smoothed state vector, C^, is smaller than the 
variance of the filtered state vector, Ck- The quadratic mean of the smoothed residual is 
closer to the measurement error (detector resolution) than the filtered one. 



4.2 Non-linear Model 

A particle's motion in a detector with magnetic field is a nonlinear process. In case of a 
non-linear system, we have to replace the transport, F^, and projection, if^, matrices in 
(fT7|) and ([IS]), respectively, by exact non-linear functions: 



Jacobian matrices of these functions (Jacobians in the following) 

djfk) d{hk) 
<9(xfc_i)' 

will be used in equations for covariance matrix propagation (fT8|) and (|20|) instead of 
and Hk, respectively. In practice, estimation with Kalman filter for a non- linear system 
shows properties similar to those of maximum-likelihood estimation: 

• The estimator is asymptotically unbiased, i.e. its bias tends to zero as the number 
of measurements increases. 

• The distribution of deviations of estimated parameters from true values approaches 
a Gaussian distribution also asymptotically, i.e for sufficiently large number of mea- 
surements. 



dm-') 



5 Particle Motion in a Static Magnetic Field 

The equation of motion of a particle with momentum p (velocity v) and charge Q in a 
static magnetic field B is: 

^ = K-Q -vx B, (26) 
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where coordinates x, y, z are in cm, p is in GeV/c, the magnetic field B is in kGauss, 
and parameter k is equah 

K = 0.000299792458 {GeV / c) kG'^ cm-\ 

The distance along the trajectory of a particle (path length) is given by: 

s = \v\ ■ t. 

The unitary vector n pointing along the direction of the trajectory is: 

dx 
a s 

Equation ( l26l) can be rewritten as: 



(27) 



dn Q ^ jpt 

— — = K- — ■nxB = K- q- nxB, 

ds \p\ 



{2t 



where q = Q/\p\- The latter equation combined with Eq. fl27|) gives a system of linear 
differential equations: 



J/' 



(29) 



dx / ds 










dy / ds 


= Uy, 








dz / ds 


= n^, 








dUx / ds 


= ^^z ■ 


Uy 


-Uy 


n 


duy / ds 


= ^^x ■ 


nz 


- ^Z 


n 


driz / ds 


= Uy- 


nx 


- ^X 


■ n 



const, 



where uji[s) = k ■ q ■ Bi{x{s)). 

6 Multiple Scattering and Energy Loss 

The ZEUS inner tracking detectors were designed using minimal material. We take ac- 
count of the effect of multiple scattering in the approximation of thin scatterers. Multiple 
scattering after traversing a material of small thickness, /, results in the perturbation of 
angles and coordinates, but the effect on the latter has an additional order of smallness 
o{l) and can be neglected. The deflection of the particle momentum p due to multiple 
scattering is decomposed into deflections in two orthogonal planes. We deflne two unit 
vectors ni, 112 which in combination with n form a right-handed Cartesian system: 



rii 



Cz X n 
\el X n\ 



1 



( -T^v \ 



V 



nx 




77-2 = riixn 



J 



1 



( 



\ 



n-y ' n,^ 



V 



with rit 



ni + nl (30) 



The direction of the momentum after the scattering is: 

n' = n + 61 ■ ni + 62 ■ n2, 



(31) 
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where 6i, 62 are random variables with 

< 6*1,2 >= 0, var (6*1 2) 



cov{9,,92) = 0. 



(32) 



Here 6ms is the well-known Moliere theory expression for RMS of the deflection angle of 
a charged particle traversing a medium [15] 



9ms {t/ Xq) 



13.6 MeV 

Pep 



t/Xo [1 + 0.038 ln(t/Xo) 



(33) 



where t/Xg is the material thickness in radiation lengths, which has to account for the 
track inclination: 

t 



I ■ + {nx/n^f + {uy/n^y 
We rewrite Eq. (!3T|) for the deflection of components: 



(34) 





( Srix ] 




' -Uy / rit ^ 




' Ux-njnt^ 


5n = 




= 01 


nx / nt 


+ 62 


ny-n-J rit 




\ Sn^ j 




\ } 







(35) 



Taking into account Eqs. fl32l) . we derive: 



< n' > 

4) 



n, 



var 

var {ri!y 
var (n^) 
cov«,ra;^) 
cov«,<) 
cov (n^, n'^) 



(^y + nln^^jn 
{nl + nlnl)/n 

^ms "-i ? 



t ) 



(36) 



^L^x riyinl - l)/nt, 
-0L riyi^'z- 



An ionization energy loss is regarded as a deterministic correction to a track energy. In 
the approximation of thin scatterer, track energy, E, after the traversal of a material is: 

E' = E- {dE/dx),on ■ t, (37) 
where {dE / dx)ion is the mean rate of ionization energy loss in the material. 



7 Specifics of Kalman Filter Implementation for the 
ZEUS Inner Trackers 

Seven equations fl2U]) describe a particle motion in a magnetic field, although five param- 
eters suffice to define the trajectory at any point. A suitable track parameterization may 
depend on the detector geometry and field shape. The magnetic field in the central part 
of the ZEUS detector is directed parallel to the z-axis. For the large part of the MVD the 
field is almost homogeneous with only a small radial component (< 1% at the edge of the 
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BMVD). For the most forward parts of the CTD and FMVD the inhomogeneity is larger, 
with reduction of the axial component by 8% and increasing of the radial component up 
to 15%. The STT detector is located outside the superconducting solenoid where the 
field is inhomogeneous. We choose a different way to proceed depending on the polar 
angle, ^, of a track (tan6' = PtlVz)'- 

• we use an option with inhomogeneous field for "forward" tracks (0 < < 60°); 

• a homogeneous field model is used for "central" tracks (60° < 9 < 120°); 

• an inhomogeneous field is used also for "rear" tracks (120° < 9 < 180°). 

The set of measurements, {m^}, with its covariance matrices, {V^}, and the map of 
magnetic field, 5, are input for the track fit. To develop a mathematical framework for 
Kalman filter implementation we have to make the following steps: 

• Select a convenient parameterization of the state vector, Xk- 

• Find a solution of the prediction equations, fk{xk-i), and a function to project the 
vector Xk to the measurement, hk{xk). 

• Obtain Jacobians of latter functions 

d{fk) djhk) 
d{xk-i)' d{xk)' 

• Define covariance matrix of the process noise, Qk- 



8 Cylindrical Parameterization for central tracks 

The magnetic field at the central region of the ZEUS superconducting solenoid is nearly 
parallel to the z-axis (5^., By 0) and has almost constant strength. Therefore we 
approximate it as homogeneous on the path from one point to the next. The system of 
equation (129|) looks as 



dx / ds 




dy / ds 


= Uy, 


dz / ds 


= n^, 


drix / ds 


= iO^-Hy 


duy / ds 


= —Uz ■ n. 


driz / ds 


= 


Q 


= const. 



where Uz = n - q ■ Bz. The component is constant and the angle (azimuthal), 0, of the 
track direction with the x-axis depends linearly on s: 

(f){s) = 00 - ^^zS, 

n^{s) = rit cos{(f)o - Uzs), , , 

ny{s) = rit sin(0o - ^^zs), 

nz{s) = Uzo, 
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where (f)o,nzo are initial values at s = 0. A pair of conserved quantities can be derived 
from (|38|): 



x(s) H ny(s) 

y{s) n^{s) 



1 

Uq Hm, 



(40) 



with initial values, xq, yo, "^xo, ""-j/O- Coordinates can by expressed via the track direction: 



1 / ^ 1 

x[s) =Xq ny[s) H nyo, 

1 1 

y{s) = H nx{s) n^^o- 



(41) 



In a homogeneous field, the particle trajectory is a helix. For the case of axial (cylindrical) 
symmetry, a natural replacement of particle coordinates, x and are the radius, r, and 
the r^^-coordinate at radius r, which we denote as u. The relation between these pairs 
of parameters reads: 

u 

(42) 



X T' COS 'r^ , 

y = r sm ^, 



and 



r 



r arctan ^ = 2r arctan ■ 



2r arctan 



r — X 



(43) 

With the usage of an arc- length in the xy-plane, st, corresponding curvature u and 
parameter A = cot 6 (cotangent of the polar angle of the particle direction) 



X 



n-. 



(44) 



(45) 



St - s-rit, ^ - rit' 
we obtain the solution for particle coordinates: 

xit) = Vq cos sm(0o ~ t) -\ sm (pQ, 

y(t) = To sm 1 — cos((/)o —t) cos 00, 

Tq UJ UJ 

z{t) = Zo + —t, 

UJ 

t = W-Su 

where ro, uq are values at t = 0. The particle which is located at a radius, ro, given t = 0, 
then arrives at a radius, r, given the value of t, which satisfies the equation: 

= Tq + T + S* sin a — (5* sin a + T) cos t + S cos a sin t, 



2 



(46) 



UJ ' 



re- 



solutions of the latter equation are 



^1,2 



2 arctan 



5* cos a 

U -2T -2S sina 



D = 



Ar(2ro + Ar), Ar = r — tq. 



D ■ (D - 2T- 2g sing) 



5^ cos a 



(47) 
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The solution t2 (with minus sign) corresponds to a shorter path length. We describe 
a particle in a homogeneous magnetic field by a state vector at a reference cylindrical 
surface of radius r^: 

xl = {uk, Zk, (j)k, Afc, qk) , (48) 

where 

Uk = ry9-coordinate at radius r^, 
Zk = -2-coordinate, 

(pk = angle of xy-projection of track direction with the x-axis, 
Xk = cot 6 at radius rk, 

Qk = Q/Pk, inverse momentum signed according to particle charge, Q. 

Such cylindrical parameterization looks natural for the barrel tracking detectors. An 
analogous state vector was used for the implementation of the Kalman filter formalism 
for the ALEPH Time Projection Chamber [12j . 



8.1 Cylindrical Parameterization: Prediction Equations 



In the prediction stage of the Kalman filter, the state vector Xk is propagated at the next 
reference radius, r^^i = + Ar^. We obtain this transformation from 



where 



Uk+l 
Zk+l 

4>k+i 

Qk+l 

Xk+1 

Vk+l 
^k 



2rfc+i arctan 
Zk + —tk, 

4>k — tk, 

Qk, 



Uk I . , 
Tk COS sm( 

Tk ^k 
. Uk 1 , 

Tk sm 1 cos( 

rk uJk 



Vk+i 



Tk+l + Xk+l 



Tk+i arctan 



Vk+i 

Xk+l ' 



(49) 



1 + A? 



tk) -\ sin (f)k 

tk) COS0fc, 

UJk 



(50) 



— 1^ ■ Bzk ■ Qk ■ y -r 

and the variable, t^, is evaluated from (1471) . We approximate the Jacobian of this trans- 
formation as: 



d{xk+i) I d{xk) 



( duk+i/duk 
dzk+i/duk 
d(j)k+i/duk 


\ 



duk+i/d(j)k duk+i/d\k 

1 dzk+i/dcpk dzk+i/dXk 
1 d(f)k+i/dXk 
1 





duk+i/dqk \ 
dzk+i/dqk 
dcpk+i/dqk 


1 / 



(51) 



Elements of the Jacobian which always are very close to zero or unity, we set explicitly to 
or 1, respectively. We exploit the sparse structure of the Jacobian to reduce computations, 
as will be discussed in Sect. [TTl Nontrivial elements of the Jacobian are presented in 
appendix A. 
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8.2 Cylindrical Parameterization: Projection of State Vector 
to MVD Measurement 



The origin of the local coordinate system of a MVD sensor is given by the vector f\,- The 
unit vector, r?, is perpendicular to the sensor plane. We define the axis of measurement 
by the unit vector, m, which is located in the sensor plane and is perpendicular to strips. 
A state vector Xk is defined at a cylindrical reference surface of a radius, r^. We can 
define the radius, r^, in such a way that the reference point will be close to the sensor. In 
the immediate vicinity of the reference point, we linearize equations fl49f5UI) with respect 
to the variable, tk- 

Xitk) =Xfc + — COS0fc, 

y{tk) =i/fc + — sin0fc, (52) 

Z{tk) =Zk+l + —tk. 

UJk 

A condition of the trajectory intersection with the sensor plane reads: 

[{r{tk)-f,)-n]=0. (53) 
The variable advance, At^, to travel from the radius, r^, to the sensor plane is: 



At, = -h. 



ak =S|cos0, + ^sin0, + g|A,, (54) 

h = {xk - Xc) + {Vk - Xc) Uy + {zk - Zc) n^. 

To obtain the expected measurement, hk{xk), we project the position vector in the local 
frame, f{Atk) — rc, to the measurement axis, fh: 

hk{xk) = [{r{Atk) - fc) ■ m] 

= 4^ + i^k - Xc) + {yk - Vc) rriy + {zk - z^) m,, ^^^^ 

Ck = rrir, cos (pk + my sin (pk + m-^\k. 
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Elements of the Jacobian, d{hk)/d{xk), are: 



dhk/duk = ^ dAtk/duk - m^f^ + ruy^ 



dhk/d(j)k = ^dAtk/d(j)k + ^^i-rn^sm(j)k + myCos(j)k), (56) 

dhk/dqk = 0, 
with derivatives of Atk 

dAtk/duk =i (n.ff-n,f^ 



dAtk/d(j)k = -^jjz (^x sin0fe - n^^ cos0fc) • 



(57) 

( T) Gin /A, — ti nriG rh, i 

To exploit the sparse structure of the Jacobian and reduce computations we approximate 
the Jacobian for specific cases: 

d{hk)ld{xk) = ( 1^ 1 H it ) , for 



(58) 



d{hk)ldM = f §^ 



dhk n dh 



^ I , for ^ 0. 



8.3 Cylindrical Parameterization: Projection of State Vector 
to CTD Measurement 

Each sense stereo wire runs at a small angle, a, and its location in the xy-plane at 
coordinate z is: 

w = fy, + {z - Zc) r'yj, (59) 

where z^ is the ^-coordinate of the nominal center of the CTD . A "planar drift" approx- 
imation is used to render measurements in space [13]. Drift distance is measured along 
the "planar drift measurement axis", m: 



=-ny/\nl 
ruy = +n^/\n\, 



(60) 



which is obtained by rotating the vector, n, through +90°. The vector n depends linearly 
on the z coordinate: 

n = + {z - Zc)p'^. (61) 

A state vector Xk is defined at a cylindrical reference surface of a radius, r^. We define 
the radius, r^, in a way that the reference point is close to the point where the trajectory 
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hits the planar drift plane. Close to the reference point, we use linearized equations of 
motion (|52|) . A condition of the trajectory intersection with the planar drift plane reads: 



[{f{tk)-w)-n] = 0. (62) 

The variable advance, Atk, to travel from the radius, r^, to the planar drift plane is a 
solution (of smallest absolute value) of a quadratic equation, (At^)^ + Atk bk + Ck = 0: 



Atki,2 = i {-bk ± ^Jbl-Aak Ck) , (63) 



with coefficients 

flfc = ^kxP'wx + "^kyP'wy, 

bk = Akx Vkx + Bkx P'wx + ^ky Vky + Bky p'^y, 
Ck Bkx T^kx ~l~ Bky T^ky i 

Akx = (cos0fe - Xkr'^^)/uJk, Aky = {sincpk - Xkr'^y) /uJk, 

Bkx ■^k ^wx (^fc ^cj^wx^ Bky yk ^wy (^fc ^cj^-u 



(64) 



Vkx = [Pwx + (Zk - Zc)pIJ UJk/Xk, Vky = Pwy + (Zk - Z^)^'^y UJk/Xk- 



UIJ/' 



wy 

The expected measurement, hk{xk), is the drift distance. To evaluate it, we project the 
position vector in the planar drift system of the wire, f{Atk) — w, to the measurement 
axis m: 

hk{xk) = [ir{Atk) -w) -rn]. (65) 

To "stretch" the projected value according to the stereo angle, a, we have to replace m 
by m/ cos a in the following formulas. The expected measurement is a linear function of 
the Atk'. 

hk{xk) = Atk Ck + frixBkx + niyBky, ^gg^ 



Ck = {m^Akx + rUyAky) / iOk- 

We approximate the Jacobian, d{hk)/d{xk), by setting its elements which are very close 
to zero or unity, explicitly to or 1: 



(67) 



dihk)/dixk) = ( 1 n n it ) ' i^^^i ^ ^0"' 

d{hk)/dixk) = 1 1^ ) , for \Atk\ < IQ-^. 
Nontrivial elements of the Jacobian are defined in appendix B. 

^ ^ — * — * 

The axial wires of the CTD run parallel to the z-axis and parameters r'y^ and p'^ 
vanish in (159|) and (16T|) . respectively. A condition of the intersection of the trajectory 
with the "planar drift plane" results in Eq. [621 which has the solution 

Atk = -bk/ak, 

ak = {cos (j)k Pwx + sin (l)k Pwy)/ ^^k, (68) 

bk (-^fc ^wx) Pwx ~l~ (Z/fc ^wy) Pwy 
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A measurement vector for an axial wire, m^, is either one-dimensional (drift distance) 
or two-dimensional (drift distance and z position). Let's consider the vector of expected 
measurement, hk{xk), for a general, two-dimensional case 

with the first component (drift distance) and second (z position), which are defined in 
(|65|) and fl52l) . respectively: 

(70) 



hki (xk) = (xfc + ^ cos 0fc - r^^) m^^ -h (y^ + ^ sin 0^ - r^j^) m^y, 
hk2{Xk) = Zk + ^^tk- 

We approximate the Jacobian, d{hk) / d{xk) as: 

/ 1 ^ \ 

= dhk, 1 9/.., J' l^^^l ^ 10"^ (71) 

and 

/ 1 \ 

d{hk)ld{xk) = dhk2 1 ' l^^'^l < ^^^^ 

Elements of the Jacobian are presented in appendix B. 

8.4 Cylindrical Parameterization: Process Noise 

We evaluate the components of a vector of particle direction, n, using parameters 0, A: 

We obtain deviations of parameters 0, A, induced by multiple scattering, from Eq. (1351) : 



50 = Oi VT+A^, 5A = -^2 Vr+A2, (74) 

where ^1,^2 are random variables defined by fl5^ . Nonzero elements of the matrix, 
describing multiple scattering in one scatterer, are: 

Qh = ^nisi^ + ^^)-> QaA = (1 + A^), (75) 

with RMS of the deflection angle, Omsi which is defined by Eq. f l33|) . The matrix, Qk, in 
Eq.f llSI) takes into account a summary effect of multiple scattering: 

Qk = ^FikQiFf^, with Fik = d{xk)/d{xi), (^^6) 

i 

and, therefore the index i runs over all elements of material on the path from {k — 1)*'' 
to k^'^ state. 
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9 Cartesian Parameterization in an Inhomogeneous 
Magnetic Field 

The following choice of track parameters at a reference z-coordinate is suited for forward 
tracks (n^ > 0): 

(x, y, t^, ty, g), (77) 



X 



where 

X = x-coordinate in the Cartesian coordinate system of ZEUS, 

y = y-coordinate in the Cartesian coordinate system, 

tx = nx/uz track slope in xz-plane, 

ty = Uy/uz track slope in yz-plane, 

Q = Q/\p\j inverse momentum signed according to particle charge, Q. 

This parametrization will be called "cartesian" . The implementation of the Kalman filter 
technique in an inhomogeneous magnetic field is analogous to those described in [16] . In 
the following we discuss the case of forward tracks. The rear tracks are specified in 
Subsec. 19.71 

9.1 Cartesian Parametrization: Equations of Motion 
in Inhomogeneous Magnetic Field 

For forward tracks we can use the z coordinate as independent variable instead of the 
path length in Eqs. fl29l) . The equations rewritten w.r.t. z coordinate read: 

dx/dz = tx, 
dy/dz =ty, 

dt^/dz = q ■ K ■ Aa:{tr,,ty,B), (78) 

dty/dZ = q ■ K ■ Ay{tx,ty, B), 

q = const. 



where the functions A^^Ay are 



A, = {l + tl + tl)-^ 



ty-{txB, + Bz)-{l+tl)By], 
-t, ■ {tyBy + Bz) + (1 + tl)Bx 



(79) 



To transport track parameters in the inhomogeneous field from plane Zq to plane z, we 
solve the latter equations with initial values defined at zo 



Xo = (xo, 2/0, txo, tyo, go)- (80) 



Three methods are used to solve Eqs. (|78l) . depending on the distance, s = z — z, 
between these planes. 

1. |s| < 10 cm: a parabolic expansion of the particle trajectory is used 



Ol 
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x{z) = Xo + t^o ■ s + ^ ■ qo ■ K ■ ■ s"^ 
y{z) = yo + tyo- s + I- qo- K- Ay ■ 
tx{z)= t^Q + qo- A^- s, 



ty{z) = tyO+qO- Ay- S, 

Qiz) = go- 

2. 10 cm < |s| < 60 cm: the classical fourth-order Runge-Kutta method [I^ is 
selected to find the solution of the equations fl78l) . 

3. |s| > 60 cm: a fifth-order Runge-Kutta method with adaptive step size control 
IT4l is used. 



9.2 Cartesian Parametrization: Equations for Derivatives 

The Jacobian of transformation of parameters given at zq to z, d{x)/d{xo), is defined as: 



d{x)/d{xo) 



( 1 




V 





1 








dx 
dy 

1 

dty 





dx 
dy 

dtyO 

dtyO 
1 





dx \ 
dy 
dt^ 
dty 

1' J 



(82) 



Elements of the latter Jacobian which are very close to zero or unity, are set to or 1, 
respectively. Nontrivial elements of the Jacobian fl82l) for short distance < 10 cm) we 
approximate as: 



dx/dt^o = s, 

dy/ dt.M = ^qQKs^ 

dt^/dtyo = qQKs^ 

dx/dqo = 2 sM^., 
dt^/dqo = KsAri 



dAy 

xQ 

yo ' 



1 2 dAj, 



dx/dtyQ 
dy/dtyo = s, 

dA 

dty/dt^o = qoi^s-^ 

dy/dqo = ^Ks'^Ay, 
dty/dqo = Ks A,, 



^3) 



xO 



'-yi 



with derivatives dA^/dtyo and dAy/dt^o, which we define below. 

To swim derivatives at long distance (|s| > 10 cm), we define equations for deriva- 
tives as described in [TB] and solve them by a Runge-Kutta method simultaneously with 
equations of motion. The magnetic field is smooth enough even in the STT area and, 
therefore we regard Eqs. fl78p as almost invariant with respect to small shifts by x and y. 
Derivatives with respect to initial Xq, y^ are trivial : 

dx'^/dxo = (1, 0, 0, 0, 0), 
dx^/dyo = (0, 1, 0, 0, 0). 
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To obtain equations for dx/dt^o, we differentiate equations flTHj) with respect to t^o and 
cliange the order of the derivative operators d/dtxo and d/dz on the left hand sides : 

d/dz{dx/dtM) = dt^/dt^o, 
d/dz{dy/dtxo) = dty/dt^, 

d/dz{dtJdt,o) = qq-k- [{dAJdtx){dtx/dtxo) + {dAx/dty){dty/dtxQ)] , (84) 
d/dz{dty/dtxo) = qo- K - [{dAy/ dtx){dU/ dUo) + {dAy/dty){dty/dtxo)] , 
dq/dt^ = 0, 

where 

dAJdt, = U ■ Ax/{1 + tl + tl) + (1 + tl + tl)^ -{ty B,~2-t,- By) , 
dA^/dty = ty ■ AJ{1 + tl + tl) + (1 + tl + tl)^ ■ it, ■ + B,) , 

dAy/dt, = tx ■ Ay /{I +tl+ tl) + {l+tl+tl)-2 ■ i-ty ' By " B ,) , 
dAy/dty = ty ' Ay / (l + + t2 ) + (1+^2^^2)1. - B y + 2 ■ ty ■ B x) . 

Initial values for the solution of latter equations are: 

dx^'/dt.o = (0, 0, 1, 0, 0). (85) 

The equations for dx/dtyo are analogous to Eqs. (18^ , but the initial values are : 

dx^/dtyo = (0, 0, 0, 1, 0) . 

To obtain equations for dx/dqo, we differentiate Eqs. ( ITHI) with respect to go and change 
the order of the derivative operators d/dqo and d/dz in the left parts : 

d/dz{dx/dqo) = dt^/dqo , 

d/dz{dy/dqo) = dty/dqo , 

d/dz{dtJdqo) = K-A^ + K-qo- [{dAJdtx){dtJdqo) + idAJdty){dty/dqo)] , (86) 

d/dz{dty/dqo) = k ■ Ay + k ■ qo ■ [{dAy/ dtx){dtx/ dqo) + (dAy/dty) (dty/dqo)] , 

dq/dqo = 1 . 

Initial values for the solution of latter equations are : 

dx'^/dqo = (0, 0, 0, 0, 1). (87) 

9.3 Cartesian Parameterization: Projection of State Vector to 
MVD Measurement 

To project a state vector (1771) to a BMVD measurement we use the method described 
in Subsect.[H21 The state vector, Xk, is defined in the reference plane with coordinate, 
z = Zk- We locate the reference plane close to the MVD sensor and, therefore use a linear 
expansion of the trajectory: 

■^(^k) ■^k ~l~ ^xk ^k 

y(sk) =yk + tyk Sk (88) 

z(Sk) = Zk + Sk- 
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A condition of the trajectory intersection with the sensor plane reads: 

[(r(sfc)-f,)-n] = 0, (89) 

where fc and n are the origin of a local MVD sensor system and the unit vector which 
is perpendicular to the sensor plane, respectively. The variable advance, As^, to travel 
from the reference plane at Zk to the sensor plane is: 

Ac. = -h. 

ak =txknx + tykny + n^, (90) 
h = {xk - Xc) + {Vk - Xc) Uy + {zk - Zc) n^. 

Analogous to Eq.[55l we obtain the expected measurement, hk{xk), by projecting the 
position vector in the local frame, f{Ask) — Tc, to the measurement axis, m: 

hk{xk) = [(fiAsk) - rc) ■ rh] 



Ask 



Ck + {xk - Xc) + {yk - Uc) my + {zk - z^) 



UJk X'^K ■^cj ""x I \yK ycj ""y i x'^K '^c) ""zi (91) 

Ck ^xk ~l~ miy tyk ~\~ fTL^. 

Nontrivial elements of the Jacobian, d{hk)/d{xk), are: 



dhk/ dxk = rux - CkUr, / ak, dhk/dyk = rriy - CkUy / ak, 
dhk/dt^k = dhk/dxk ■ Ask, dhk/dtyk = dhk/dyk ■ Ask- 



(92) 



Derivatives with respect to slopes have an additional order of smallness o{Ask) and we 
approximate the Jacobian, d{hk)/d{xk), for the BMVD: 



a("o/9(io =(ii §^ It It » ) ■ i^^'i s 1° 



-3 

5 



d{hk)/d{xk) = ( 1^ ||j ) , for |A.,| < 10- 



(93) 



Sensors of the FMVD are almost perpendicular to the 2;-axis and, therefore rix^y ~ 0. 
We locate the reference plane at the position of the FMVD sensor {Ask = 0). Taking 
into account latter remarks, we obtain from Eq. (p2|) the Jacobian, d{hk)/d{xk), for the 
FMVD 

d{hk)/d{xk) = { my ) . (94) 



9.4 Cartesian Parameterization: Projection of State Vector to 
CTD Measurement 

The linear expansion of a particle trajectory fl88p defines the particle coordinates in the 
immediate vicinity of a stereo wire. An approach to obtain the projection of cartesian 
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state vector to CTD stereo measurement is similar to those discussed in Subsec. 18.31 A 
condition of the trajectory intersection with the planar drift plane reads: 



[if{sk)-w)-n]=0, (95) 

where the coordinate of the wire, w, and vector, n, are defined by fl59l) and fl6Tl) . respec- 
tively. The variable advance, As^, to travel from the reference plane to the planar drift 
plane is a solution of a quadratic equation, (As^)^ + Ask bk + Ck = 0: 

Asfc = ^ [-bk + ^Jbl-Aak Ck) , (96) 

with coefficients 

bk = Akx Vkx + Bkx V'wx + ^ky Vky + -Bfcy P^wy 
Cfc ^kx T^kx ~l~ -Sfcy T^kyi 



^kx tkx ^wx^ ^ky tky ^wy^ 

Bkx -^k ^wx (^/c ^c) "^wxi Bky Uk ^wy (^fc ^c) ^wyi 

'Pkx Pwx ~l~ (^fc ^c) Pwxi "^ky Pwy ~l~ (^^k ^c) Pwy 

We obtain the expected measurement, hk{xk), i.e. drift distance, as in fl65l) : 



(97) 



hk(xk) 



Xk + tkx Ask - r^x - {zk - Zc + Ask) r'^x 
Uk + tky Ask - r^y - {zk - Zc + Ask) r'^y 



my, 



(98) 



where the m have to by replaced by m/ cos a to take into account the stereo angle, a. 
Derivatives with respect to slopes have an additional order of smallness o{Ask) and we 
approximate the Jacobian, d{hk)/d{xk), for the CTD stereo measurement: 



dhk dhk dhk dhk 

tyk 



d{hk)/d{xk) = ( g S ) , for |A.,| > 10- 



d{hk)/d{xk) = ( 1^ 1^ ) , for \Ask\ < W 



(99) 



Elements of the Jacobian, d{hk)/d{xk), are presented in appendix C. 

Axial wires of the CTD run parallel to the 2;-axis and parameters, r'^^ and p'^ vanish 
in (l59l) and (I6T]) . respectively. A condition of the intersection of the trajectory with the 
planar drift plane leads to Eq. (p5|l . which has the solution: 

Asfc =-bk/ak, 

^k txkPwx ~l~ tykPwyj (100) 

bk i^Xk 1"wx) Pwx ~l~ (Z/fc ^wy) Pwy 

We consider the vector of expected measurement, hk{xk), for the general, two-dimensional 
case 
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with the first and second component being a drift distance and z-position, respectively: 

) m.^ + {Vk + tky 



hkli^k) — {Xk + tkx ^Sk 

hk2ixk) = Zk + Ask. 



We approximate the Jacobian, d{hk)/d{xk), as: 

dhki dhkA 



d{hk)/d{xk) 




dhki 
(^Vk 
dhk2 
<^yk 



and 



d{hk)/d{xk) 



( dhki 

dXk 

dhk2 
\ dxk 



dh 



k2 



dhki 



'dTkx 
dhk2 

dtkx 








dtky 

dhk2 

dtky 








for I As J > 10 



-3 



for lAsfcl < 10" 



(102) 



(103) 



(104) 



where we take into account an additional order of smallness o{Ask) for derivatives with 
respect to track slopes. We obtain nontrivial elements of latter Jacobians in appendix C. 



9.5 Cartesian Parameterization: Projection of State Vector to 
STT Measurement 

Signal wires of a given STT layer are arranged in a plane perpendicular to the 2;-axis 
with coordinate z = z^. We locate the reference plane at the position of the layer, i.e. 
Zk = The particle trajectory inside a straw tube we approximate by a straight line. 
The latter line and the signal wire are described as lines which pass through points fk 
and fw and have directions Hk and respectively: 



(105) 



Components of the vector of particle direction, n^, we calculate using track slopes tkx, tky- 
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rikx 



tkx 



'-kx ''ky V ~^ ky V ^fca; "r '^ky 

The expected measurement is a drift distanc^ in the straw, which is evaluated as a 
distance between these two lines: 



1 + tl^. + tf.y 



hk{xk) 



[rk - ry,) -nkXn^ 
\nk X n^l 



(107) 



^ We expect that left-right ambiguity of the drift distance is resolved and, therefore regard it as a 
signed value. 
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After simple calculations the expected measurement reads: 

hk{xk) = ^ ^ — • (108) 

y 1 ~l~ (j^kx ^wy ^ky '^wx) 

The Jacobian of the latter transformation we can approximate as: 

5(MMx.) = (|| II O), (109) 

with 

dh}^l 'IT'wy I \/l ~l~ ^kx f^wy ^ky ^wx) j 

(110) 

dhk/dyk = n^x / \/l + (tkx n^y - tky n^xf. 

9.6 Cartesian Parameterization: Process Noise 

We evaluate deviations of track slopes induced by multiple scattering from Eq. ( |35l) : 

Stkx 



yukzJ rikzrikt rif^^nkt 

where 61,62 are random variables defined by (l32l) . Nonzero elements of the matrix de- 
scribing multiple scattering in one scatterer are: 

Qtxtx = ^ms (1 + tlx) (1 + ^Ix + tky)y 

Qtyt,=6La+tly){l+tl+tl^), (112) 

Qta: ty = 6^g tkx 4y (1 + tkx ^'^ky)^ 

with RMS of the deflection angle, which is defined by Eq. f l5^ . The matrix, Qk, in 
prediction equation ( ITSl) has to account for a summary effect of multiple scattering on a 
path from {k — 1)^^ to k^^ state, and is therefore evaluated analogous to (!76|) . 

9.7 Cartesian Parameterization for Rear Tracks 

For rear tracks {n~, < 0) we use a parameterization analogous to those for forward tracks. 
The meaning of parameters x, y, q is identical with fl77|l . For rear tracks we define slopes 
w.r.t. negative direction of the 2-axis: 

tx = -nx/n,, ^^^g^ 

ty ^y I '^z • 

Equations of particle motion for rear tracks are identical to flTSl) for forward tracks, but 
with slightly different definition of functions A^iAy-. 



Ax = {l + ti + tl) 

Ay = {l+tl+tl) 



ty{-txBx + B,) + {l + tl)By\, 

-tx-{-tyBy + B,)-{l+tl)Bx]. ^ ^ 

Equations 0881) for linear and (ISTj) for parabolic expansions of trajectory can be used for 
rear tracks also, if we regard the expansion w.r.t. z-coordinate decrement, s = zq — z. 
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10 Global Parameterization 



A global perigee parameterization of tracks [13] is used for analyses in the ZEUS experi- 
ment. The perigee parameters are parameters of a helix, which are defined at the track's 
point of closest approach to the 2;-axis: 

= Q/Rh, QDh, zh, cote) , (115) 

where 

(pH = angle of xy-projection of track direction with the x-axis, 

Q/Rh = helix curvature signed by a particle charge, Q, 

QDfj = signed minimal distance to 2;-axis, 

zh = ^-coordinate at point of closest approach, 

cot^ = cotangent of track direction w.r.t. 2;-axis. 

Transformations between local parameters (cylindrical or cartesian) and global ZEUS 
perigee parameters are given in appendix D. 

11 Fast Computations with Kalman Filter Technique 

Most of the calculation by the Kalman filter technique is in the following procedures: 

• transportation and projection of track parameters fl2^ and evaluation of Jacobian 
matrices (125|) : 

• matrix operations in prediction fllSp . filter fl2ip and smoother fl23p equations; 

• search of a track crossing with material to evaluate effects of multiple scattering 
and energy loss. 

Approaches to fast computation with Kalman filter technique were discussed for the 
magnet tracking [I7j,[18] at the HERA-B detector. 

To reduce computations we use a flexible strategy for propagating track parameters 
and derivatives in the inhomogeneous field, as described for forward and rear tracks in 
Subsec lom For long (s > 10 cm) distances we use numerical integration of the equations 
of motion, but integrate derivatives together with a "zero trajectory" that allows to 
reduce computations. For short distances (s < 10 cm) we use parabolic expansion (|8T|) 
of the particle trajectory, which is very fast in computations. 

To keep the computational effort at a minimum we exploit the sparse structures of 
the Jacobian matrices. The Jacobian of track propagation includes elements which are 
very close to or 1, therefore we use Jacobian approximations and set such elements to 
or 1. The Jacobians for cylindrical (!5T]) and cartesian (!82|) parameterization contain 
only 11 and 10 nontrivial elements, respectively. To calculate the product of matrices 
Fk Cfc_i -FJ in fllSp we implement functions, which take into account a sparse structure 
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of the matrix Fk- For example, the function for 10 nontrivial elements of the implies 
73 multiplications, which is much smaller than 200 multiplications needed for the case of 
the completely filled matrix F^ of size 5 by 5. 

The Jacobians of projection transformation, if^, are approximated also, as shown in 
(158|) . (1671) . (I7T]) . (!72l) etc. We implement corresponding functions for the calculation of 
products of C^-^Hl and (1 - KkHk)Ct^ in ([H]) or (1 - HkKk)Vk in These 
functions take into account the sparse structure of the matrix if^,. For example, only 20 
multiplications are sufficient to obtain the matrix {1—Kk Hk) Cj^~^ for the option with one 
nontrivial element in the matrix Hk- This has to be compared with 100 multiplications 
needed for the completely filled matrix of size 5 by 1. 

To evaluate the effects of multiple scattering and energy loss, we describe the distribu- 
tion of material in the ZEUS inner trackers by using about 1800 separate volumes. After 
crossing a given volume, a particle can reach only a limited number of other volumes. We 
implement an approach called volume navigation [19J for fast search of a track's crossings 
with these volumes. Using the Monte Carlo technique, we evaluate for each volume a 
list of volumes, which can be crossed subsequently. On average, one list includes about 
7 subsequent volumes. The lists are used to navigate a fast search of track crossing with 
volumes. 

The described approaches have been programmed ^20j in C++. We follow recipes of 
effective programming of numerical calculations [H] and implement STL containers to 
store objects like hits, states, tracks etc. 

Table 1: Computing time of the track fit per ZEUS event on a PC with processor Intel 
CPU 3.06GHz for different groups of tracks. 



Fitted tracks 


Fraction 


Field model 


Time / event 


Forward {6 < 60°) 


59% 


inhomogeneous 


12 ms 


Central (60° <9 < 120°) 


23% 


homogeneous 


7 ms 


Rear (9 > 120°) 


18% 


inhomogeneous 


1 ms 


All tracks in event 


100% 


(in) homogeneous 


20 ms 



A ZEUS event contains up to 100 fitted tracks and about 30 tracks on average. The 
longest tracks include about 80 hits in the central area, 50 hits in a transition region 
and 30 hits in the very forward direction. Fitting all the tracks in one event takes 20 ms 
on PC with processor Intel CPU 3.06GHz (see Tabled]) and 46 ms with processor Intel 
CPU IGHz. About of 77% of tracks are fitted using the inhomogeneous field as shown in 
Tabledl The computing time for these tracks is comparable with those which are fitted 
using the homogeneous field approximation. 
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Figure 1: Standard deviations of the pull distributions of fitted track parameters in perigee 
parameterization (jllSp for MC simulated muons versus the momentum. 
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The precision of fitted parameters depends on the resolution and details of the per- 
formance of the ZEUS trackers and will be discussed in the next note [20]. Here we 
would like to mention that approximations implemented to reduce computations, are not 
made at the expense of track parameter precision. Evaluation of the covariance matrices 
of fitted parameters in (12T!) and (123!) are the most complicated computations, including 
operation with the transport, projection and process noise matrices. Calculated vari- 
ances of these matrices are in good agreement with the residuals of the fitted parameters. 
Standard deviations of pull distribution (residuals normalized by their estimated error) 
are close to unity for different track momenta, as shown in Fig.[U 

12 Conclusions 

We consider a mathematical framework for the rigorous approach to a common track fit 
using the trackers in the inner region of the ZEUS detector: CTD, BMVD, FMVD and 
STT. We discuss track models and likelihood functions in such a multi-component tracker. 
The approach offers a rigorous treatment of field inhomogeneity, multiple scattering and 
energy loss. The track fitting procedure makes use of the Kalman filter technique. 

We describe details of the mathematics for the fast implementation of a Kalman filter 
for the cylindrical drift chamber, barrel and forward silicon strip detectors and straw drift 
chambers. The cases of homogeneous and inhomogeneous field are considered. 

We discuss how to reduce computations and make the track fitting procedure fast. 
Average computing time of track fitting in one ZEUS event is about of 20 ms on a PC 
with processor Intel CPU 3.06GHz. 
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13 Appendix A: Jacobian of prediction transforma- 
tion in cylindrical parameterization 



We use 



derivatives of tk to calculate the Jacobian flSTl) of prediction transformation fH9l) 



dtk/duk = 



COS(0fc -—) - COs(0fc 

rk 



dtk/dzk = 0, 



Tk COs(0fc Ifc, 



ifc) H sm tfc 

rk 



dtk 

dtk/dX, 



k — 



dtk/duk = 



= -rk dtk/duk, 

dtk/duk, 

^k 

2(1 — cos tk) + Tk oJk sin(0 ) — sin( 

L Tk 



Nontrivial elements of the Jacobian are: 



9 / \ 

Tk cUfc cos(0 tk) +ujk sm tk 

rk 



Sk : 



duk+i/duk = 

^ rk+i 
Vk+i 



duk+i/dXk 



cos ^ + ^ sin(0fc - tk) dtk/duk 



sin^- ^ cos(0fc - tk) dtk/duk 



duk+i/d(j)k = Jj!'xk+i [si^'^'fe ~ sin(0fc - tk) (1 - dtk/d(j)k)] 

= ul'y'k+i [~ + cos(0fe - tk) (1 - dtk/d(t)k)] 



rk+i 

^k Xk+1 



(116) 



1 + A 



^ (cos(0fc - tk) - cos0fc) + sin(0fc - tk) dtk/dXk 



rk+i 

^k Vk+l 



-r4\^ (sin(0fc - tk) - sin0fc) - cos(0fc - tk) dtk/dXk 



duk+i/dqk 



dzk+i/duk 
dzk+i/dXk 

d(t>k+i/duk 

d(j)k+i/dqk 



[cos^fc - cos(0fc - tk) + ujk sin(0fc - tk) dtk/duk] 



Qk 



rk+i 



Qk ^k Vk+l 

^^dtk/duk, 
tk 

-dtk/duk, 
-dtk/dqk- 



[soKpk - sin(0fc - tk) - uJk cos(0fc - 4) dtk/duk] , 



dzk+i/d(j)k= ^dtk/d(l)k, 
XkUJk 



dzk+i/dqk 
d(f)k+i/dXk 



^-dtk/duk, 

i + Aj. 



(117) 
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14 Appendix B: Jacobian of projection transforma- 
tion for the CTD in cylindrical parameterization 

Elements of the Jacobian fl67|) of projection transformation f l66|) for the stereo CTD are: 

dhk/dzk 
dhk/d(f)k 

dhk/dXk 



Ck dAtk/ dzk - r'^^ - ruy r'^y, 
Atk 



Ck dAtk/d<Pk 



iOk 



rrix sin0fc + niy cos (j)k) 



Ck dAtk/dXk - Atk 



XkCk 
T+Al 



+ 



ujk (j^^ ~'~ ^'wy) 



with derivatives 

dAtk/dzk 
dAtk/d(j)k 

dAtk/dXk 
dbk/dzk 
dck/dzk 
dak/d(f)k 
dbk/d(f)k 
dak/dXk 
dbk/dXk 

dck/dXk 

dAkx/dXk 

dAky/dXk 



Atk dbk/ dzk + dck/ dzk 
2 Atk dk + bk ' 

Atfc dbk/dcpk 

Atfc2 dak/ d(j)k + 2 Atk ak + bk 

Atk dbk/dXk + dck/dXk 
Atfe2 dak/ dXk + 2 Atk dk + bk' 

i^-^kx Pwx ~^ -^ky Pwyj Pwx ^wx Pwy ^wyi 
{Skx Pwx "I" -^ky Pwy) ~ T^kx ^wx ~ '^ky f'-u 

^ (- sin (Pk p'wx + cos 0fc v'w)) , 
^ (-sin (pk Vkx + cos (pk Vky) , 
p'^^ dAkx/dXk + p'wy dAky/dXk, 
Vkx dAkx/dXk + Vky dAky/dXk 

^ {Akx Vkx + Aky Vky) , 



1 


Xk 




Xk{ 


1 + AD 


1 


Afc 


^k 


.1 + A^ 


1 


Afc 




1 + A^ 



rricos (pk 



~T~ J^ky I ky ) 1 

^kr'wx) + r'wx 



'^k^wy) "I" ^wy 



Elements of the corresponding Jacobian (ITTIl for the axial CTD look as: 

(-A4sin</), + ^ cos0fc) 

Q A _ 

Ik 



dhki/dcpk 



dhk2/duk 
dhk2/d(j)k 
dhk2/dXk 
dAtk/dcpk 



^k 
X 



+ ^(AtfcCos0fc + ^ sm(/)fcj 



ak 



k (r, "Mk — n ^k \ 



Xk dAtk 
^k d(j)k ' 

Atk 



i^k 
At, 



ak i^k ^^"^"^ 



:ii8i 



(119) 



(120) 



Sm(/)fc -Pwy COS (pk) 
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15 Appendix C: Jacobian of projection transforma- 
tion for the CTD in cartesian parameterization 

Nontrivial elements of the Jacobian fl99|) of projection transformation fl98|) for the stereo 
CTD are: 



dhu/dxk = m, + M^ dhk/dyk = rriy + 

dhk/dt^k = 1^ ■ Asfc, dh/dtyk = ^ ■ Ask, 

with 

2 Ask ak + bk 2 Ask au + bk 

= m^^ (tkx - r'^x) + ^wy {tky - r'^y)- 
Elements of the corresponding Jacobian (llU3p for the axial CTD read as: 



dtyk 




dhk2 


_ dAsk 






dhk2 


_ dAsk 


dtyk 


(^Vk 



with 



dAsk/dxk = dAsk/dyk = 

ak ttk 

M^yj fTT'wx ^kx ~l~ ^^wy ^ky 



(121) 



(122) 



16 Appendix D: Conversions from Local to Global 
Parameters 



Track parameters, Uq, Zo,(j)o, Xo,qQ, at the beginning of a central track, which is fitted 
using the cylindrical parameterization (HHl) . are converted to perigee parameters (1115^ : 



Q/Rh = KB,qo^l + Xl 
QDh - 



where 



cote 



tH 



- To sin(— - 00 + tn) + (cost// - 1) / (Q/Rh), 

zo + XotH / (Q/Rh), 
Xq, 



(123) 



arctan 



I /{Q/Rh) - rp sin(Mo/ro - 0o) 
ro cos(Mo/ro - 0o) 



f sign(g/i?^f). 
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We convert fitted cartesian parameters at the beginning of a forward track ,xo, yo, tox, toy, Qo, 
into perigee parameters flllSp : 

(pH = arctan 



where 



Q/ Rh = i^B^qo- 



Q_ 



2 I +2 ' 



-1 + 3^0 sin 0j^ + A'o cos 0/^ 
= QAR^ ' ^^^^^ 

zh = Zo + 



tlx + tov Ql Rh 



cote 



H 



tOX + t^y 

toy 



$0 = arctan 

tox 

Xo =-yoQ/RH+ °^ 



tox + tn 



''Ox T I'Oy 

yo = xoQ/Rh+ ^ ■ 



The transformation from cartesian to perigee parameterization for rear tracks is similar 
to fll24p for parameters (pn, Q/Rh, QDh, but differs for parameters zh and cot 9h- 

Zh = zq- 



cote 



tox + tlyQ/Rn (125) 



H 



tox + toy 
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